Deep Inelastic Pion Electroproduction at Threshold 
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We consider the cross section of the deep inelastic pion electroproduction on a proton target at 
threshold for in the region 5 — 10 GeV^. The corresponding amplitudes are described in terms 
of two form factors which we calculate using light cone sum rules (LCSR) to leading order in QCD 
and including higher twist corrections. Our results suggest a considerable change from a small 
region that can be treated in the soft pion limit using current algebra. In particular, we obtain a 
tt'' to 7r+ producton ratio of order 1/3 and significant nucleon helicity-flip contributions. 
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I. INTRODUCTION 

In recent years there has been increasing atten- 
tion to hard exclusive processes involving emission 
of soft pions in the final state. One reason is that 
such processes often provide the main background to 
the simpler reactions that one wants to study, and 
one would like to estimate the background as precise 
as possible. The main motivation is, however, differ- 
ent. It turns out that hard production of soft pions is 
interesting in its own right and can provide one with 
new insights in the hadron (nucleon) structure. The 
novel physical feature is in this case the presence of 
three disparate scales: the large momentum transfer 
Q, the QCD scale A = Aqcd, and the pion mass 
771^ which goes to zero in the chiral limit, so that 
Q ^ AqcD ^ TOir- Moreover, the limits — 
and Q ^ oo do not commute. If, parametrically, 
<^ A'^/m^r then the standard approach based on 
chiral symmetry and current algebra is applicable. 
In this way the celebrated soft pion theorems arise 
[l[ : the amplitudes with pion emission arc calculated 
in terms of the amplitudes without pions induced by 
a chirally-rotated current, and pion emission from 
the external hadron lines. There exists vast litera- 
ture on this topic, see e.g. 0, d, 3, extending the 
original KroU-Rudcrman prediction to include cor- 
rections in the pion mass and/or the pion (photon) 
momentum. 

In the opposite limit ^ A^/itIt^ these tech- 
niques do not apply, and the accepted general phi- 
losophy is to use QCD factorization to separate con- 
tributions of large momenta of the order of the hard 
scale Q in pcrturbatively calculable "coefficient func- 
tions" in front of nonperturbative matrix elements 
that contain contributions of smaller scales. The 
methods of current algebra (or chiral perturbation 
theory) can then be used to simplify the matrix el- 
ements rather than the full physical amplitude in 
question. The present study is fuelled by the ob- 



servation made in Ref. Q that in the asymptotic 
Q — > oo limit the pion electroproduction can be de- 
scribed in the framework of the usual pQCD fac- 
torization formalism for exclusive processes [1, 0, Q 
and involve an overlap integral of chirally rotated 
nucleon light-cone distribution amplitudes (DAs). 
The particular combination of the DAs involved in 
the pion electroproduction turns out to be differ- 
ent from that in the electromagnetic (or weak) elas- 
tic nucleon form factors, so that different compo- 
nents of DAs potentially could be separated from 
the comparison. This is interesting, since nucleon 
light-cone DAs remain to be poorly known. Simi- 
lar ideas have been applied in a somewhat different 
context in Refs. [ol. IToj. 

The essential requirement for the applicability of 
the pQCD factorization is a high virtuality of the 
exchanged gluons and also of the quarks inside the 
short distance subprocess. The main problem is a 
numerical suppression of each hard gluon exchange 
by the Us/tt factor which is a standard perturba- 
tion theory penalty for each extra loop. If, say, 
ag/iT '--^ 0.1, the pQCD contribution to baryon form 
factors is suppressed by a factor of 100 compared to 
the purely soft term. As the result, the onset of the 
perturbative regime is postponed to very large mo- 
mentum transfers since the factorizable pQCD con- 
tribution has to win over (nonfactorizable) nonper- 
turbative effects that are suppressed by extra powers 
of but do not involve small coefficients. There 
is a growing consensus that such "soft" contribu- 
tions play the dominant role at present energies. In- 
deed, it is known for a long time that the use of 
QCD-motivated models for the wave functions al- 
lows one to obtain, without much effort, soft con- 
tributions comparable in size to experimentally ob- 
served values (see, e.g. [llllil, A current fash- 
ion [ij, [Hi is to use the concept of generalized par- 
ton distributions (GPDs) to describe/parametrize 
soft contributions in various exclusive reactions, see 
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[H, [13, [ill for recent reviews, and the models of 
GPDs usually are chosen such that the experimen- 
tal data on form factors are described by the soft 
contributions alone, cf. Rcfs. [H, [13, |2l|. The 
soft pion elcctroproduction was considered in this 
framework in [l^, [11]. A subtle point for these 
semi-phenomenological approaches is to avoid dou- 
ble counting of hard rescattering contributions "hid- 
den" in the model-dependent hadron wave functions 
or GPD parametrizations. 

Another approach to the calculation of form fac- 
tors for moderately large is based on the light- 
cone sum rules (LCSR) [HI, [11| . This technique is 
attractive because in LCSRs "soft" contributions to 
the form factors are calculated in terms of the same 
DAs that enter the pQCD calculation and there is 
no double counting. Thus, the LCSRs provide one 
with the most direct relation of the hadron form fac- 
tors and distribution amplitudes that is available at 
present, with no other nonperturbative parameters. 

In this paper, we suggest to use light-cone sum 
rules to calculate the amplitudes of pion elcctropro- 
duction. The basic object of the LCSR approach is 
the correlation function 

j dx e-^-^- (iV(P)7r(fc) |r{j (x)ry(0)}|0) 

in which j represents the electromagnetic probe and 
is a suitable operator with nucleon quantum num- 
bers. The final state nucleon and the pion are ex- 
plicitly represented by the state vector (A^(P)7r(fc)|, 
see a schematic representation in Fig. [T] When 




FIG. 1: Schematic structure of the light-cone sum rule 
for pion electroproduction 

both the momentum transfer and the momentum 
[P'Y = (P — q + k)^ flowing in the ry vertex are large 
and negative, the asymptotics of the correlation 
function is governed by the light-cone kinematics 
^ and can be studied using the operator prod- 
uct expansion (OPE) T{j{x)T]{0)} ~ J2C^{x)O^{0) 
on the light-cone x^ = 0. The ^^-singularity of 
a particular perturbatively calculable short-distance 



factor Ci{x) is determined by the twist of the rele- 
vant composite operator Oi, whose matrix element 
(A^(P)7r(fc)|Oi(0)|0) is given by an appropriate mo- 
ment of the (complex conjugated) pion-nucleon gen- 
eralized DA. Using current algebra and chiral sym- 
metry these matrix element can be reduced to the 
usual nucleon DAs. Next, one can represent the an- 
swer in form of the dispersion integral in (P')^ and 
define the nucleon contribution by the cutoff in the 
invariant mass of the three-quark system, the so- 
called interval of duality sq (or continuum thresh- 
old). The main role of the interval of duality is that 
it does not allow large momenta > sq to flow 
through the 77- vertex; to the lowest order 0{a^) one 
obtains a purely soft contribution to the form fac- 
tor as a sum of terms ordered by twist of the rele- 
vant operators and hence including both the leading- 
and the higher- twist nucleon DAs. Note that, in dif- 
ference to the hard mechanism, the contribution of 
higher-twist DAs is only suppressed by powers of 
|(P')2| ~ 1-2 GcV^ (which is translated to the 
suppression by powers of the Borel parameter after 
applying the usual QCD sum rule machinery), but 
not by powers of Q^. This feature is in agreement 
with the common wisdom that soft contributions are 
not constrained to small transverse separations. 

The LCSR expansion also contains terms generat- 
ing the asymptotic pQCD contributions. They ap- 
pear at proper order in as, i.e., in the 0(as) term for 
the pion form factor, at the order for the nucleon 
form factors, etc. In the pion case, it was explicitly 
demonstrated [H, [13 that the contribution of hard 
rescattering is correctly reproduced in the LCSR ap- 
proach as a part of the 0{as) correction. It should 
be noted that the diagrams of LCSR that contain 
the "hard" pQCD contributions also possess "soft" 
parts, i.e., one should perform a separation of "hard" 
and "soft" terms inside each diagram. As a result, 
the distinction between "hard" and "soft" contri- 
butions appears to be scale- and scheme-dependent 
[26j . During the last years there have been numer- 
ous apphcations of LCSRs to mesons, see [H, [1^ 
for a review. Following the work [sl] nucleon elec- 
tromagnetic form factors were further considered in 
this framework in Refs. [Sll, [H, [H, [s^ and the weak 
decays At p^vt, Ac Mvi etc. in [sl, [H, [37| . 
A generalization to the NjA transition form factor 
was worked out in (38| . 

The presentation is organized as follows. Section 2 
is mainly introductory. We introduce here the neces- 
sary notations, explain the kinematics and define the 
two independent invariant amplitudes (generalized 
form factors) that are relevant for the pion electro- 
production at threshold. The main nonperturbative 
input in our calculations is provided by the general- 
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ized pion-nucleon distribution amplitudes which are 
calculated using current algebra in Section 3. The 
light-cone sum rules are constructed in Section 4. 
The numerical analysis is done in Section 5 which 
contains our main results. Finally, Section 6 con- 
tains a short summary and conclusions. 



II. KINEMATICS AND NOTATIONS 



the same frame, the electromagnetic current is de- 
fined as 



j'in^) = euu{x)jf_i,u{x) + edd{x)j^d{x) , (2.5) 

and P{W) is the kinematic factor related to the 
c.m.s. momentum of the subproccss 7* ((7)-)-p(P') 
7r(fc) + N{P) in the final state: 



In this paper we consider pion electroproduction 
from a proton target 

e(/') + p{P')^e{l) + n+{k) + n{P) , 

e(O+p(P'He(0 + 7r0(fc)+p(P). (2.1) 

Basic kinematic variables are 

q = l'-l, s={l' + P'f, W^^{k + Pf, 
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" — 



P' -q _W^ + Q'^-m% 



P' ■ I' 



(2.2) 



"N 



The identification of the momenta is clear from 
Eq. (j2.ip : toat is the nucleon and 

ttitt the pion mass, 
respectively. In what follows wc neglect the elec- 
tron mass and the difference of proton and neutron 
masses. 

The differential cross section is equal to 



da 



1 



\M\' 



dydQ'^d(t) (27r)5 64(s - m%) 
Here M is the amplitude: 



^f3{W)d<j).^d{cose) . 

(2.3) 



M = -^u{lMl'){Nn\f;-mp) , (2.4) 

6 and (j)^^ arc the polar and azimuthal angles of the 
pion in the final nuclcon-pion cm. frame, respec- 
tively, (f) is the azimuth of the scattered electron in 
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(mAf-l-TO^) 



/ 4 V" W'^ 



1^2 



m) 



w 



(2.6) 



Alternatively, one could use the Mandelstam 
t-variable (of the 7*p — > TriV subprocess) 
t = {P' -Pf: 



dt = 2|A:i||A:/|d(cos6') , 



(2.7) 



where ki is the c.m.s. momentum in the initial 
state: 



M/2 



(2.8) 

We will be interested in the pion production ex- 
actly at threshold, in which case the pion momen- 
tum is simply proportional to that of the final state 
nucleon: 

k^ = 5P^, (5 = TO^/TOAr~0.15, (2.9) 

and assume that the photon virtuality is large, 
^ ^Qcr>- The conditions of Lorentz covariance 
and electromagnetic current conservation imply that 
the amplitudes at threshold can be parametrized in 
terms of two generalized form factors which we de- 
fine as 



(iV(P)^(fc)b— (O)lp(P')) = -f A^(P)75 |(7m9' - 1, i) ^Gl'^m - !^GJ^(Q2)1 n{P') . (2.10) 

I 



Hereafter N{P) is the usual Dirac u-spinor, = (DIS) cross section from unpolarizcd protons close 
93 MeV. For example, the S'-wavc contribution to to the pion-nucleon threshold is given in terms of 
the structure functions of the total deep-inelastic the generalized form factors as: 
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f3{W) + (2mjv + TO^)^ 



(47rA) 



m%{mN + m^) 



(2.11) 



The calculation ol Gl^ {Q^) and GJ^(g2) in the 
light-cone sum rule approach presents our main goal. 

Having in mind the practical construction of light- 
cone sum rules that involve nucleon DAs, we define 
a light-like vector by the condition 



q • z = 0, 







(2.12) 



and introduce the second light-like vector 

V,=P.-\-,^.^ P' = 0^ (2-13) 



so that P ^ p a the nucleon mass can be neglected, 
toat 0. The photon momentum can be written as 



P-q 
P ■ z 



(2.14) 



We also need the projector onto the directions or- 
thogonal to -p and z, 



pz 



and use the notation 



dp = a^ip^ 



(2.15) 



(2.16) 



for arbitrary Lorentz vectors a^. In turn. a± denotes 
the generic component of orthogonal to z and p, 
in particular 



p-q 

q±^L = q^i . 

p ■ z 



(2.17) 



We use the standard Bjorken-Drell convention 
[sot for the metric and the Dirac matrices; in par- 
ticular, 75 = i7°7^7^7^, and the Levi-Civita tensor 
e/ifAcr is defined as the totally antisymmetric tensor 
with eoi23 = 1- 

Assume for a moment that the nucleon moves 
in the positive direction, then p"*" and z~ are 
the only nonvanishing components of p and z, re- 
spectively. The infinite momentum frame can be 
visualized as the limit p+ ^ Q — > oo with fixed 
P ■ z = p ■ z ^ \ where Q is the large scale in the 
process. Expanding the matrix element in powers of 
l/p+ introduces the power counting in Q. In this 



language, twist counts the suppression in powers of 
Similarly, the nucleon spinor N{P, A) has to be 
decomposed in "large" and "small" components as 



N{p,xy 



1 



(^^+^Ti)N{P,X) 



2p- z 

--N+{P,X) + N-'{P,\), (2.18) 



where we have introduced two projection operators 



A+ 



2p- z 



A" 



2p-z 



(2.19) 



that project onto the "plus" and "minus" compo- 
nents of the spinor. Note the useful relations 

liN{P)=mNN+iP), /N{P) ^'^^N-(P) 

rriN 

(2.20) 

that are a consequence of the Dirac equation 
fN{P) = mArA^(P). Using the explicit expressions 
for N{P) it is easy to see that A+iV = 7V+ ~ ^/p+ 
while A~N = N- ^ 

Note that all expressions are invariant under the 



reparametrization 



azn where a is a real num- 



ber; we will use this freedom to set equal to the 
"minus" component of the distance between the cur- 
rents in the operator product. 



III. PION NUCLEON GENERALIZED 
DISTRIBUTION AMPLITUDES 



A. Leading Twist 

Our purpose will be to get a complete classifica- 
tion of three-quark generalized distribution ampli- 
tudes (GDAs) of a riTT"*" and p7r° pair in the limit 
that the momentum of the pion relative to that of 
the nucleon is small which will be a generalization of 
po| . This can be done using soft-pion theorems, in 
the spirit of the work Q . In order to develop the nec- 
essary formalism we first consider the leading twist 
GDA in some detail. 
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1. Nucleon Distribution Amplitudes 

Wc begin with quoting the necessary portion of 
the definitions and results from ,40]. 

The notion of hadron distribution ampUtudes 



J 



(DAs) in general refers to hadron-to- vacuum matrix 
elements of nonlocal operators build of quark and 
gluon fields at light-like separations. In this paper 
we will deal with the three-quark matrix element 



(0| e^^''u]^{aiz) [aiz, aoz]^,^^ (a2z) [a2Z, aoz]j, j d'^ (032) [032, aoz],,, j. A)) 



(3.1) 



r 



where \p{P, A)) denotes the proton state with mo- 
mentum P, — M'^ and helicity A. u, d are 
the quark-field operators. The Greek letters a,/?, 7 
stand for Dirac indices, the Latin letters A: refer 
to color, z is an arbitrary light-like vector, ~ 0, 
the ai are real numbers. The gauge-factors [x, y\ are 
defined as 



, y] = P exp 



ig j dt{x-y)f,A''{tx+{l-t)y) 




(3.2) 

and render the matrix element in (|3.ip gauge- 
invariant. To simplify the notation we will not write 
the gauge-factors explicitly in what follows but im- 
ply that they are always present. 

Taking into account Lorentz covariance, spin and 
parity of the nucleon, the most general decomposi- 
tion of the matrix element in Eq. ()3.1|1 involves 24 
invariant functions To the leading twist-three 

accuracy only three of them are relevant. We can 
write in a shorthand notation 

4 (0| e^-'"'=<(aiz)^.^^(a2z)d^(a32) b(P, X)),^,,,^, = 

(3.3) 

where 

iv^U,,=i^C)^p (75iV+)^ 

{tiU,,={^'Jj^pC)^f, (7^75^+)^ (3.4) 

stand for the Lorentz structures. 

In turn, , A^, Tf can be written as 

F{a^p-z) = J Vxe-'P'^' '''''' F{xi), (3.5) 

where the functions F{xi) depend on the dimen- 
sionless variables Xi, < Xi < = 1 which 

correspond to the longitudinal momentum fractions 



carried by the quarks inside the nucleon. The inte- 
gration measure is defined as 

Jvx = J d.Tida;2dx3 5(xi -I- a;2 + — 1) ■ (3.6) 



Note that 



(wi)Q/3,7 = (wi)^a,7 , 
(ai)a/3,7=— (ai)/3a,7 7 
{tl)al3,i = {tl)pa,j ■ 



(3.7) 



Since the operator on the l.h.s. of Eq. (|3.3p is sym- 
metric under the exchange of the two w-quarks, this 
property implies that the V- and T-functions are 
symmetric and the A-function is antisymmetric un- 
der the exchange of the first two arguments, respec- 
tively: 

T/f(l,2,3)=Ff(2,l,3), 
A?'(l,2,3)=-Af(2,l,3), 
rP(l,2,3)=rP(2,l,3). (3.8) 



In addition, the matrix element in Eq. (|3.3p has to 
fulfill the symmetry relation 



<(1)«^^(2)<(3)|P) 



<(1)«^,(3)4(2) \P) 
+ {Q\e^^'^u\{3y^{2)d';,{l)\P)=Q (3.9) 

that follows from the condition that the nucleon 
state has isospin 1/2: 

+ 1)) (0| £'^'=<(1)^.^^(2)4(3) |P) = , 

(3.10) 

where 



t4 



(3.11) 
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and T± are the usual isospin step-up and step-down 
operators. Applying the set of Ficrz transforma- 
tions 

{vi)',p,a=^ {vi - ai - ii)„_g ,y 

(<l)7/3,a = - (wi +ai)„^_^ (3.12) 

one ends up with the condition 

2rf (1, 2, 3) = [Vr - AP] (1, 3, 2) + [Vr - A?] (2, 3, 1) , 

(3.13) 

which allows to express the tensor DA of the leading 
twist in terms of the vector and axial vector distri- 
butions. Since the latter have different symmetry, 
they can be combined together to define the single 
independent leading twist-3 proton DA 

mXi,X2,X3) = [V[-AP]{xi,X2,X3) (3.14) 



which is well known and received a lot of atten- 
tion in the literature. The neutron leading twist DA 
$3 (xi, a;2, aJa) can readily be obtained by the inter- 
change of u and d quarks in the defining Eq. p.3p . 
For all invariant functions F = A,T proton and 
neutron DAs differ by an overall sign: 

FP(1,2,3) = -F"(1,2,3), (3.15) 

as follows from the isospin symmetry. This property 
is retained for all twists. 



2. Pion-Nucleon Generalized Distribution Amplitudes 
of twist-3 

Our aim is to describe the generalized distribution 
amplitudes (GDA) of a pion-nucleon system with 
small invariant mass in the similar formalism. To 
this end we define for the pTr^-system: 



4{0\e''''ul,iaiz)uj,{a2z)d';{a3z) |p(^^, A)^°)^^_3 = (75)7^7^ \vr\vi)c.,3.s + Af {m)^p^s + Tf" {h)o^i3,s 

(3.16) 

and, similar, for rnr^: 

4 (0| e'^''ul{a,z)u^{a2z)d';ia3z) \n{P, X)7t+)^^_^ = (75)7*7^ kr^(«i)a/3,5 + ^^^(01)^/3,5 + T^^hU^s 

Jir L 

(3.17) 

An extra 75 is needed to conserve parity. Similar to the proton case, the symmetry of the two w-quarks 
implies that the GDAs V and T are symmetric, and A is antisymmetric to the exchange of the first two 
arguments, respectively. 

As observed in j^, the GDAs can be calculated using current algebra for small invariant masses of the 
pion-nucleon system in the soft pion limit. One obtains 



^> A',/' 

(3.18) 

^' A',/' 

(3T9i 



(0|O|7r'^(fc)iV/(P, A)) = -j^mQt 0]\Nf{P, A)) - E (0|O|iV^'(P, A')) N{P, A') j^^^r'^f. fN{P, A) 

and similar for the pion-nucleon final state 

(^°(fc)^/(^,A)|Ot|0) = -l.(^N}{P,\)\[Qh0^m + ^^Jp^ J2 ^(-P'^) lil^rf^'N{P,X'){Nj,iP,\')\O^\0). 

A' ,/' 



r 



Here = 93 MeV is the pion decay constant defined axial charge 
as (0|g7M75 5T>|7r''(A:)) = iS^'^f^k,, and gA ^ 1.25 is 
the axial charge of the nuclcon. O is the nonlocal 
three-quark operator and Qg is the operator of the 



d^a;g(a;)7o75yg(a:^) 



(3.20) 
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where the t° are the usual Pauh matrices. The sec- 
ond term in Eq. (|3.18[) . Eq. (|3.19p corresponds to 
the pion bremsstrahlung from the ingoing (outgo- 
ing) nucleon. This term corresponds to a ttN sys- 
tem in P-wave and can. in principle, be separated 
from the first (commutator) contribution that is S- 
wave by considering the angular distributions in the 
TriV system. Note also that the bremsstrahlung con- 
tribution vanishes at the threshold of pion produc- 
tion Wth = -A/ + but becomes significant for 
W — Wth ^ "m-n- This contribution is determined 
in terms of the nucleon DAs that eventually com- 
bine to produce the nucleon electromagnetic form 



factors; it can always be added. In what follows we 
will concentrate on the S-wave (commutator) term. 

For our purposes we need to specify O to be the 
relevant three-local light-ray operator. 

O ^ Olf^ = £'^''=<(aiz)w^^(a2z)d^(a3^) (3.21) 

Here the superscript (uud) indicates the flavor 
content and simultaneously the order in which 
the quark fields are positioned on the light ray. 
For example, O"^^ and O^^" correspond to the 
Ua{aiz)ui3{a2z)d^{a^z) and Ua{aiz)di3{a2z)u^(a^z) 
configurations, respectively. One obtains 



{Q\OX\n{PA)^^{k))^^{{0\Ol%\n{PA)^\k^^ 

=--^[{Q\[Ql Olf;\ |n(P, A)) + *(0| [Qi , Olf^]^?, A)) } + . . . 

(o|oj:;^^,b(p,A)7r"(fc))=-f (o|[g^,o:;5',]b(^,A)) + ... (3.22) 

where the ellipses stand for bremsstrahlung contributions. Calculation of the commutators [Q5, 0^^;^(2/)] 
and [Q\,0'^Ji{y)] reduces to a chiral rotation 



[Qh<lf] = -[\] 75<?/'. (3.23) 

For the three-quark operators one uses the chain rule [A, BCD] = [A,B]CD + B[A,C]D + BC[A,D\ to 
obtain: 

[01^0^%'^^]= ^{(75)qaOa!^^ + (75)/3aOa7q - (75)7aO"^a} ) 
1 



!5 



.Olt,]^--{{^,UOlf, + (75)^a0^5;' - il^hxOTpi] ■ (3.24) 



Taking the nuclcon-to- vacuum matrix clement, Eqs. (|3.3p . (|3.16p and (j3.17p . and using symmetry relations 
p.7|) and Fierz transformations (|3.12|) one obtains after some algebra 



AT\i,2, 3)=i={iFi"(l, 3, 2) - iyi"(2, 3, 1) - ^^^(2, 3, 1) -f ^^^(l, 3, 2) + T^{2, 3, 1) - T^{1, 3, 2)}, 
(1, 2, 3)=-i={iv^i"(l, 3, 2) + ^^^(1, 3, 2) + ri"(l, 3, 2) + \v^{2, 3, 1) + \a^,{2, 3, 1) -f T{\2, 3, 1)}, 
Tr^(l,2,3)=^{AlH2,3,l) + AJ'(l,3,2)-Fi"(2,3,l)-T/i"(l,3,2)}, (3.25) 



and 



Fr°(l,2,3)=iF/'(l,2,3), 
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(l,2,3)=-^?(l,2,3), 

Tr"(l,2,3H^rf(l,2,3), (3.26) 
which is the desired result. One can further simplify (|3.25p using the substitution in Eq. (|3.13p to obtain: 
Vr^il, 2, 3)=i={Fi"(l, 3, 2) + Vni, 2, 3) + Fi"(2, 3, 1) + A^l, 3, 2) + AU2, 3, 1)}, 

(1, 2, 3)=-^[vr{3, 2, 1) - Vni, 3, 2) + A'/(2, 1, 3) + (2, 3, 1) + A^S, 1, 2)}, 

Tr\l, 2, 3)=^{a»(2, 3, 1) + 3, 2) - ¥^{2, 3, 1) - 3, 2)}. (3.27) 

I 



Note that the functions V, A, T for the both cases 
riTT"*" and pi:^ satisfy the symmetry conditions (j3.8p 
which follow directly from their definitions in (|3.16p , 
p.l6[) . respectively. On the other hand, the isospin 
relation similar to (j3.13p is not valid, since the pion- 
nucleon pair can have both isospin 1/2 and 3/2. 

B. Higher Twists 

Taking into account Lorentz covariance, spin 
and parity of the nucleon and the pion, the most 



general decomposition of the three-quark matrix 
element in Eq. (|3.ip involves 24 invariant func- 
tions. In total, there are three leading- twist- 

3 invariant functions (Vi,Ai,Ti), nine of twist- 

4 (F2,A2,y3,^3,5i,Pi,T2,T3,r7), nine of twist-5 
( ^4 , ^4 , , ^5 , £"2 , ^^2 , T4 , T5 , Tg ) , and three of twist- 
6 (Vfi, ^6, Tfj). Using the shorthand notation of 
Ref. [40| for the relevant Lorentz structures, we de- 
fine 



4- (0|£'^'=<(aiz)M^^(a2z)4(a3z)|iV(P,A)^(fc)) = 



(75)7*-^ 



+ n"'^(«5)a/3,5 + V^'^{vq)o,P^S + AJ^(ai)a/3,5 + ^j'^(a2)a/3,5 

+ 2^3^("3)"^'.'5 + 2^4^('^4)q/3,<5 + ^5^(05)0/3,15 + Aq^ {aQ)al3,S 



+Tl itl)al3,S + T2 (,t2)af3,S + T^ (^3 )q;3,(5 + {t4)al3,6 



(3.28) 



The calculation of higher-twist pion-nucleon GDAs much more cumbersome. The complete set of the 
is similar to the leading twist, but proves to be necessary Fierz identities for the Lorentz structures 



9 



{vi)ai3,s, {ai)ai3,s {U)a0,s IS given in Ref. [40|. The 
results for the GDAs of all twists are collected in 
Appendix A below. In addition, in Appendix B wc 
present the results for the functions V^^ and Ai^ 
that appear in the oflF- light-cone 0{x^) corrections 
in the OPE involving three-quark currents, 

4(0|e'J'^<(aix)u^^(a2x)d^(a3x) |7V7r) = 



(75)757^1 



Vi + 



Ti 



(ai)a/3,5 



{vi)al3,S 



ttN 



(3.29) 



for the case that the light-cone positions of two of the 
three quarks coincide, to the first subleading order 
in the conformal expansion. The function T/^ docs 
not contribute to the sum rules to our accuracy. 



IV. LIGHT CONE SUM RULES 

A. The soft pion limit 

For technical reasons, it is convenient to write 
the sum rules for the complex conjugated amplitude 
with the pion-nucleon pair in the initial state. To 
this end we consider the correlation function 



T;^(F,g) 



(4.1) 



{0\T{^,{0)fr{x)}\N{P)7T{k)), 



where 



(Ohp(O)liV(P)) =A?m^A^(P), (4.2) 

is the so-called loffe proton interpolating current 
jilj l , Af is the corresponding coupHng. 

In the limit |fc| for fixed and (F')^ = 
{P + k — q)^ the correlation function in (14. ip can 
be calculated using PCAC and current algebra in 
terms of the correlation functions without a pion 



and involving chirally-rotated currents 

r;^(p,g) = 

=-j- I J d^xe^'^^0\T{[Ql,VpmJr{^)}\N{P)) 



I / d^xe^'^^Q\T{r^p{Q)[Ql,jl^{x)]}\N{P)) 



-\- bremsstrahlung terms , (4-3) 
where Qg is the axial charge and the bremsstrahlung 
contributions correspond to pion absorption by the 
initial state nuclcon. 

The commutators can easily be evaluated with the 
result 



QifJ^{x)]^Q 



(4.4) 



where r± - ^(r^ ± zr^), Q± = ^{Ql± ^Ql) etc. 
Also 



[Qt ,Vp{x)]=---j=lBrin{x) 
[Q|,77p(x-)]=-^75ryp(a;) , 



(4.5) 



where 77„ is the neutron current 

r^n{x) = -e*^" [d\x)C^^d^{x)\ -i^^uHx) . (4.6) 
The contributions of interest to (|4.ip are those sin- 
q q 



a) 



p 



p 



k ^. 



b) 



FIG. 2: Schematic structure of the pole terms in the 
correlation function (14.111 



gular in the vicinity of P'^ m^, see Fig. [2j Note 
that in addition to the nucleon pole, Fig. [2^, one 
has to take into account the semidisconnected con- 
tribution with the pion-nucleon intermediate state. 
Fig. ^Bp. For example, for n7r+ one obtains 



T; "(P,g) = 



'TV 



pi2 



75 { {luq^ - qy 



viJ.q^ Q-mr^ 



'AT 



2ra 



N 



N{P) 
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-75 



- (P' - kf 



(4.7) 



where the elhpscs stand for less singular contributions. and F2 arc Dirac and PauU electromagnetic 
neutron form factors, respectively 



{N{p'm^mN{p)) = N{p') 



N{P) 



(4.8) 



and A„7r+ is the coupling of the loffe current to the mr^ state: 



1 7?p(0) \N{P' - fc)7r(fc)) = -_A,wmw75iV(P' - k) . 

Jtt 



(4.9) 



On the other hand, using the representation in Eq. (|4.3p instead, one obtains 



rp7: + n(p . i\\mN 1 mN+ f" 



2mN 



GpiQ')-^^GT{Q')}l5N{P) 



2'mN 



FHNiP) + 



(4.10) 



where the form factors in the first line are defined as 



{N{P')\A,{0)\N{P)) = 7V(F') 7.G^(Q') - 7^Gp{Q') - i:^Gt{Q') l5N{P) . (4.11) 



2mN 



2 m 



Note that Gt{Q^) ~ because of the isospin and 
CP invariance. 

In the same approximation the pion-nucleon cou- 
pling is given by A„^+ = (-73) K (cf- Eq. g31)) 
so that the terms in the second line in Eqs. (|4.7p 
and (|4.10p coincide: The contribution of the chiral 
rotation of the nucleon current is identically equal 
to the semidisconnected contribution in Fig. ^Bp of 
the pion coupling to the current. Equating the re- 
maining contributions in the first line of Eqs. (j4.7[) 
and (|4.10p one obtains 

2 

Q^GT* {Q^) = ':^Ga{Q^) + 0{m^/K, \k\/K) , 

GT^ (O')=0 + 0(m,/A, |fc|/A) , (4.12) 

which is the classical result [H, [3j S 3 ■ The consider- 
ation of the pii^ state is similar, the only difference 
being that the commutator of the electromagnetic 
current with the axial charge vanishes in this case, 
so that the both amplitudes vanish at threshold, up 
to corrections in m^r and/or |fc|: 

Q^Gf{Q^)=Q + 0{m,IK\k\/K), 

"(02)^0 ^ 0(m,/A, |fc|/A) . (4.13) 

To avoid misunderstanding, note that our 



(schematic) derivation does not include 
bremsstrahlung contributions, which have to 
be added. The derivation breaks down for large 
momentum transfers, however. One way to see 
this is that the limit fc — > that is implied in the 
current algebra relations, becomes far away from 
the physical region. In particular, the position of 
the pole in the pion-nucleon intermediate state 
in Fig. [2)d (P' — kY — mj^ moves away from 
the nucleon pole at (P')^ = mj^. Requiring that 
(P' - /s)2 - (P')2 < A^pp ^ m% one obtains a 
restriction for the applicability of the soft pion limit 
^ A^/m^. In what follows we develop a tech- 
nique to calculate the amplitudes in the opposite 
limit of large momentum transfers, 3> A^/rn^. 



B. The sum rules 

The general philosophy of the LCSR approach is 
explained in the Introduction. The main idea is to 
separate contributions of small distances to the cor- 
relation function in (|4.ip in terms of calculable coef- 
ficient functions in front of operator matrix elements 
and apply the soft pion techniques for the evaluation 
of the latter rather than the correlation function it- 
self. Information about the generalized form factors 
can then be extracted by matching the QCD calcu- 
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lation at moderate negative values P'^ ^ — l GeV^ 
with the dispersion representation (j4.7p in terms of 
hadronic states. 

The correlation function in (j4.ip involves several 
invariant functions that can be separated by the ap- 
propriate Lorentz projections. The structures that 
are most useful for writing the LCSRs are usually 
those containing the maximum power of the large 
momentum p+ ~ pz. We define 

A+T, = +^{pz + kz)j5 {mNA+ 4^B] N+{P) . 

JTT 

(4.14) 

The invariant functions A and B can be calcu- 
lated for sufficiently large Euclidean and P'^ = 
{P + k — q)^ in terms of the generalized distribution 
amplitudes of the pion-nucleon system using the op- 
erator product expansion. We use the following no- 
tations: 

F{x3)=J dx'^ J dxiF{xi,l- xi- x'^,x'-i) , 
F{xs) = dx'^ dx. 



X / dxi F{xi,l — xi — x^,x^) (4.15) 
Jo 

and 

F{X2)= I dx'2 I dxiF{xi^x'2T^ — Xi — x'2) ^ 

F{x2) = / dx2 i dx2 

X / dxiF{xi,X2,l ~ xi — X ), (4.16) 
Jo 

where F ~ A,V,T is a generic pion-nucleon GDA 
that depends on the three valence quark momentum 
fractions. 

To the three-level accuracy A^^ and B"^^ are 
given by the same expressions [s^l as the LCSRs 
for the proton electromagnetic form factors Fi and 
i^2, respectively, with the substitution of the nucleon 
DAs by the pion-nucleon ones: 



A 



2ed / dx3 



<3^ + 93t7^jv , 2:3 



1I 



% Jo 



4 



2e,„ / dx2{^ r dx, [^2V,-^ + 3^3^^^ + ] (^0 " ^^Vf"'''^*") 
92 Jo 



4 



I2 f>7rAr 



92 



4 '^123 



Q +92 2^N a;;m 

d ^123 



2 '"'TV 



91 



9| 



•^1345 ~ ^ "^43 ~r ^3^ 



^34 



2x2m^^/^ 

A ^ 123456 ( ' 
92 



B 



ttN_ 



2ed / dx3 



9i Jo 



V 



4 



''123 ~ ''43 



N 



2e„ / dx2 



I2 Jo 



dxi[Vi+A,r'' {x,)+ 4 

92 



N 



X2m% 

4 



''134.5 + ''193 r ^19.T — ^''43 + ^3/ 



134 



(4.17) 



where we used shorthand notations for the combina- V12345 = 2Vi — V2 — V3 — V4 — , 

tions of the DAs: 

^^123456 = -Vi + V2 + V3 + V4 + V5-Ve, (4.18) 
Vi23 ^Vi-V2-V3, A3i =A3~Ai, 

V1345 = -2^1 + V3 + V4 + 2V5 , A123 = -Ai +A2-A3, 
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^1345 = -2Ai - A3 - ^4 + 2A5 , 

^12345 = 2Ai - ^2 + A3 + ^4 - ^5 , 

Ai23456 = Ai - A2 + A3 + A4 - A5 + (4.19) 

and also 

T137 = Ti —T^ — Tj , 

J 13478 — ^Jl^ J3^J4 — J7 — -tSi 

^134678 = Ti - T3 - T4 + Te - T7 - Tg • (4.20) 

Explicit expressions for the functions y^^^^^f"^'') g^^jj 
^7r7V,A/(ji,ii) ^^^g given in Appendix B. The super- 
scripts (u) or (d) stand for the cases that one of 

I 



the M-quarks has the same position as the d-quark, 
and the two u-quarks have the same position, re- 
spectively. In addition, we use a compact notation 

q3 = q-X3iP + k), q2=q-X2{P + k). (4.21) 

In the both expressions in (|4.17p the functions with a 
"tilde" and a "hat" have X3 and X2 as an argument, 
respectively, cf. (|4.15p . (|4.16|) . Also, in the terms 
involving two integrations over the momentum frac- 
tions, the remaining momentum fraction is replaced 
by using xi -t- X2 -I- X3 = 1. 

The Borel transformation and the continuum sub- 
traction are performed by using the following sub- 
stitution rules: 



dx 



g{x) 



dx 



{q-x{P + k)f 
q{x) 



{q-x{P + k)Y Jo x2 (.s-P'2)2 



^ da; g{x) 
~ {s-P'^) 

^ dx q{x) 



1 da; 



q{x) cxp 



1 /"'da; 



xQ^ 




xAP 


M2 ) 






xQ^ 


xm\i 


xlVP 


Af2 



(4.22) 



q{xo) e" 



where M is the Borel parameter, s = ^-^pQ^ + (1 — x)m\i and is the solution of the corresponding 
quadratic equation for s = sq: 



(g2 + So - m2,)2 + 4m2^Q2 _ (q2 ^ _ 



(4.23) 



The contribution - e--""/^^' in Eq. ([422)) corre- 
sponds to the "surface terms" arising from the par- 
tial integration to reduce the power in the denom- 
inator (g — xPY = [s — P'^)'*(— a;)"* to the usual 
dispersion representation with the denominator ^ 
(s — P'^). Without continuum subtraction, i.e. in 
the limit sq ~* 00 this term vanishes. 

On the other hand, the nucleon contribution of 
interest corresponds to a pole term in the variable 
P'2 at P'2 = TO^. For the relevant projections we 
get 



A 



TV 



m% - [P' - k) 



'AT 



pi2 



'AT 



(P' - kf 



(4.24) 



where Pi and P2 arc the Dirac and Pauli proton 
(for p'k'^) and neutron (for ni:'^) electromagnetic 
form factors, respectively. The pion-nuclcon cou- 
pling A^rAf = ^ttn{P''^) is defined in Eq. (|4.9p and 



in general depends on the invariant mass P'^ of the 
pion-nucleon system. Close to threshold in the P'2- 
channel, P^ = to^(1 + 5)^, one obtains in the soft- 
pion limit A„jr+ = ~Af/\/2, Xp^^o — — A^/2. 

Note that the Borel transformation amounts in 
this case to the substitution 



1 



TO 



TV 



pl'2 



(4.25) 



'TV 



(P' - fc)2 



(1 + J)e-[™«(i+'5)'+'5Q']/A/' 



One observes that the contribution of the pion- 
nucleon intermediate state is suppressed compared 
to the nucleon contribution of interest by an ex- 
tra factor exp{(5[m^(2 + 5) + Q'^]/M'^}. Requir- 
ing that this suppression is at least as strong as 
exp{— [sq — ?ti^]/M2} in which case the pion-nucleon 
state has to be considered as a part of the contin- 
uum, we obtain 



<5[to^(2- 



S) + Q^]>so-m 



TV : 



(4.26) 
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Q^>iso-m%)/S-{2 + S)m% (4.27) 

where from > 7.3 GeV^. For smaller momen- 
tum transfers the pion-nuclcon contribution has to 
be calculated explicitly and subtracted. The prob- 



lem is that the corresponding couplings Xj^n have to 
be taken at the invariant mass of the pion-nucleon 
pair P'^ = m^(l + 6)'^ + SQ^ and this dependence 
is not known a priory. 



Neglecting the pion-nucleon intermediate state, we obtain the sum rules: 







/" 




1 

"AT 


\f 

J xa 



d.T - 



exp 



xQ2 


xm'j^ 


xAP 


A'P 






xQ^ 




xNP 





Q2 -I- x^m 



N 



(4.28) 



where the functions ^2 4(^) given in terms of the gcnerahzcd pion-nuclcon distribution amphtudes as 



X 



g^^{x)=2eA V^l^ + x I dx.V^^^ (x,) \ + 2eJ x I dxi [^2Vi + 3V3 + ^3]""^ (a;,) - V,l'^ + A^^, 



10 ) I Jo 

gl{x)^2e,^ Q%1^ + x^m%V:A + 2ejQ' (v,l^ + l^^^ 



43 + ^34 



2Trn2 f v'^^'^^^") -UT/ 
ZXrU;^ l/j -t- V 123456 



gl{x)=-2ed^j^ dxiyi"^(x0| +2e„|^ dxi[Vi + Aif'' {x,)^ , 



g\{x)=-2eaml{vl 



TN,M{d) 



V123 - V4 



43 



2eum 



N 



nN,M{u) .ttN,M{u) 



^1345 + 1^23 + ^123 ^ 2143 + A 



134 



(4.29) 



Note that the sum rules have exactly the same structure as the sum rules for electromagnetic nucleon form 
factors in Ref. [l^l, with the only difference that nucleon DAs are replaced by the pion-nucleon ones. 

For completeness, we quote also the sum rule for the axial form factor of the proton from Ref. [U (for 
the charged current defined in (|4.4p ): 



I ^ + ^2 ^,^2 ) 



xQ2 xml^X £>|(xo)e-('^"-™«)/^^' 



tAP AP 



+ XQm'j^ 



(4.30) 



with 



02ix) = { +xj dxiVaix,) j + ^x J dxi [2Ai + 3^3 + ^3] (x^) - A123 + V123 



^1345 — 2^34 -I- V43 



glix) = ^ QV123 + x-'mijV^sj + i^Q^ ( A123 + V123 ) + x'm'j^ 

2xm% ( A^^'"'' -Ai23456 



(4.31) 



r 



Note that in difference to Eq. (|4.28p this sum rule involves nucleon DAs on the r.h.s., not the pion- 
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nucleoli ones. 



RESULTS 



The sum rules in Eqs. (|On|) involve nu- 

cleon and pion-nucleon DA that can be expanded 
in contributions of conformal partial waves. To each 
order in the conformal expansion, one or several new 
nonperturbative parameters appear that can be re- 
lated to matrix elements of local operators, as de- 
tailed in Ref. [40l |. In the simplest approximation 
(asymptotic DAs) the sum rules depend on a single 
parameter which is the ratio of the twist-4 and the 
twist-3 matrix elements of the three-quark operators 
of the lowest dimension (e.g. without derivatives) 
that differ by the quark spin projection on the light- 
cone. This ratio is known to ca. 20-30% accuracy 
from the QCD sum rule calculations (see e.g. [33|). 
To the next-to-leading-order in the conformal spin, 
several more parameters enter that characterize the 
momentum fraction carried by a particular quark in 
the valence component of the nuclcon wave function, 
both in the leading and higher twists. These pa- 
rameters have been estimated using QCD sum rules 
[13, S, Hi but the accuracy of such estimates is 
suspected to be low. By this reason we prefer to use 
the set of parameters that are intermediate between 
the asymptotic and the QCD sum rule motivated 
DAs, and are chosen in such a way that the nucleon 
electromagnetic and axial form factors are described 
well within the light-cone sum rule approach, see 
Ref. 01 . We wiU refer to this set of the DAs as 
the Braun-Lenz-Wittmann (BLW) model. For lead- 
ing twist, these parameters are close to the model of 
Ref. m. 

In order to minimize further the uncertainty due 
to the choice of the DAs, the Borel parameter 
and the continuum threshold sq, we present be- 
low the results for the ratios of the sum rules for 
the pion elcctroproduction amplitudes at threshold. 
Eqs. (|4.28p . to the nucleon axial form factor, Eq. 
KM : 



and 



N 



(5.1) 



(5.2) 



For small momentum transfers, the first ratio is pre- 
dicted to be one for the ir'^n and zero for Tr^n final 
state (|4.12p . while the second ratio vanishes in both 
cases in the soft-pion limit (|4.13p . 

The LCSR (|4.30p for the axial form factor itself is 
studied in the range = i _ iq QeV^ in Ref. [H. 




. 4 
0.3 
0.2 
0.1 







V2|GJ"|/Ga 
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FIG. 3; Pion electroproduction amplitudes at threshold, 
(|2.10p normalized to the nucleon axial form factor: ir'^n 
(solid blue curves) and (dashed red curves). The 
calculation is done using the BLW model [s^l for the 
nucleon DAs. The color identification refers to the on- 
line version. 



The results are consistent (to 20% accuracy) with 
the dipole parametrization 



Ga(Q')= 



9A 



(5,3) 



with gA = 1.267 ± 0.004 and the mass parameter 
Ma = 1 GeV that is suggested by both the neutrino 
scattering and pion electroproduction data at low 

The results are shown in Fig.[3]for the range of mo- 
mentum transfers = 5 — 10 GeV^. This range is 
chosen because for lower momenta the contribution 
of the pion-nucleon intermediate state becomes large 
and must be subtracted, and, on the other hand, for 
large momenta the (uncalculated) radiative correc- 
tions 0{as{M)) are expected to play an increasing 
role. 

It attracts attention that the -Ri ratio is almost 
constant in the whole range, i.e. the Q'^GJ^ 
amplitude turns out to be proportional to the ax- 
ial form factor of the nucleon. This proportionality 
holds for all models of the nucleon DAs that we are 
considering, whereas the value of i?i does depend on 
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the shape of the DAs. We obtam 



i?^ "=0.80 ±0.10, 



i?i ''=0.45 ±0.08, 



(5.4) 



where the given error estimate corresponds to the 
difference in the sum rule predictions using BLW 
and asymptotic DAs. 

The second amplitude, GJ^, comes out similar in 
magnitude but different sign for the 7r+n and ir'^p 



final states. The sum rule results for the i?2 
shown in Fig. [3] can be parametrized as 



l + b{Q'-Ql) 



ratio. 



(5.5) 



with the reference scale Qq ~ 7.5 GcV^, the slope 
b ~ 0.1 GeV~^ and the values 



R2 "(Qo)^ 0.21 ±0.15, 
i?f p(q2)^_0.19±0.07. 



(5.6) 



where the error estimates are obtained as in 
Eq. (lOl) . 

Using our central values for the amplitudes 
we obtain the contribution to the structure function 
F2 (|2.1ip which in the sum of all terms is almost the 
same as the corresponding contribution in the soft 
pion limit, cf. Eqs. (|4.12p . (|4.13p . This is in disagree- 
ment to Ref. [H] where a much larger structure func- 
tion was obtained. The inspection shows that our 
result for the ir'^n production coincides with that of 
Q within the errors, while our prediction for the n^p 
cross section is smaller by an order of magnitude. 
As the result, the DIS cross section obtained in Q 
is completely dominated (90%) by the contribution 
of the n^p state (compare Eqs. (11) and (12) in [^), 
whereas in our calculation the tt*^ to tt"*" production 
ratio is only of order 1/3. Although the precise num- 
bers are sensitive to the models of the nucleon DAs, 
there seems to be a qualitative difference that would 
be interesting to check experimentally. 



VI. SUMMARY AND CONCLUSIONS 

We have given an analysis of the pion electropro- 
duction amplitudes at threshold in the light-cone 



sum rule approach, using the methods of current al- 
gebra to construct the complete set of pion-nucleon 
distribution amplitudes in the soft-pion limit. The 
sum rule derived in this work are tree-level and may 
be affected by large (but calculable) radiative cor- 
rections. Having this in mind, the numerical results 
presented in Section 5 have to be considered as semi- 
quantitative. Our calculation suggests that the tt^ 
to 7r+ production ratio becomes roughly of order 1 /3 
at large momentum transfers > 7 GeV^, and the 

amplitudes Gj ^ and GJ^" are not small. They are 
predicted to be similar in magnitude and have op- 
posite sign which can be tested in experiment using 
polarization transfer techniques. Working out the 
concrete experimental predictions (e.g. correspond- 
ing asymmetries) goes beyond the tasks of this paper 
and will be considered separately. 

Main assumption underlying the present calcu- 
lation is that the matching of the QCD operator 
expansion with the dispersion relation in terms of 
hadronic states can be done at sufficiently large 
scales where pions do not present independent de- 
grees of freedom. In practice, the relevant values of 
the Borel parameter are close to the cutoff scale of 
the chiral perturbation theory, so this assumption 
may need a further check and better justification. 
Another direction of further work is obviously the 
inclusion of radiative corrections to the sum rules, 
which are enhanced at large momentum transfers. 
For practical applications one may need to add con- 
tributions of pole terms and generalize our analysis 
of the pion-nucleon distribution amplitudes to small 
off-threshold pion momenta by using the Omnes rep- 
resentation. 
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APPENDICES 

APPENDIX A: PION-NUCLEON GDAS OF HIGHER TWISTS 

In this Appendix we give without derivation the expressions for three-quark pion-nucleon generahzed 
distribution amphtudes of higher twist. 

1. nir^ distribution amplitudes 

Definition of the GDAs foUows Eq. (|3.28|1 . For the twist-4 we obtain: 

(1, 2, 3)=-i=|l ({^^(1, 3, 2) - Pr(l, 3, 2) + 5^(1, 3, 2) - r3"(l, 3, 2) 

-T7"(l, 3, 2) - Kr(l, 3, 2)} + {l ^ 2}) I , 
Ar\l, 2, 3)=-^| ^ ({^3 (1. 3, 2) + P^il, 3, 2) - 3, 2) + r3"(l, 3, 2) 

+T7"(l,3,2)-T/3"(l,3,2)}-{l^2})|, 
(1, 2, 3)=-i= I i ({a^(1, 3, 2) - Pi"(l, 3, 2) + Sni, 3, 2) + r3"(l, 3, 2) 

+T7"(1, 3, 2) + V^{1, 3, 2)} + {1 ^ 2}) I , 
(1, 2, 3)==-i= I i ({ ~ A^(l, 3, 2) - Pi"(l, 3, 2) + 5^(1, 3, 2) + T3"(l, 3, 2) 

+r7"(l,3,2)-^"(l,3,2)}-{l^2})|, 
(1, 2, 3)— i ({^^(1,3, 2) + A^(l, 3, 2) + Pi"(l, 3, 2) + 3, 2) 



+2T2"(1, 3, 2) + r3"(l, 3, 2) - (1, 3, 2) - ¥,"{1, 3, 2) 
+y3"(l,3,2)}-{l^2})|, 

(1, 2, 3)=-i=|l ({ - A^{1, 3, 2) - A^(l, 3, 2) + Pi"(l, 3, 2) + 5f (1, 3, 2) 
+2T2"(1, 3, 2) + r3"(l, 3, 2) - r7"(l, 3, 2) + V^il, 3, 2) 
-y3"(l,3,2)}-{l^2})|, 
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(1, 2, 3)— ^1 ^ ({ 3, 2) + 5^(1, 3, 2) ~ T3"(l, 3, 2) + T^\l, 3, 2)} 
+ {1^2})}, 

T3"-^(l, 2, 3)-^| i ({a'2^(1, 3, 2) - A-{1, 3, 2) - Pi"(l, 3, 2) - 5^(1, 3, 2) 
+2T2"(1, 3, 2) - Tg'Xl, 3, 2) + (1, 3, 2) - ^^"(l, 3, 2) 
-V^3"(l,3,2)} + {l^2})|, 

(1, 2, 3)=-^ I i ({a^(1, 3, 2) - A^(l, 3, 2) + P^^(l, 3, 2) + 5^(1, 3, 2) 
-2T2"(1, 3, 2) + T^{1, 3, 2) - T^{1, 3, 2) - V^{1, 3, 2) 

-F3"(l,3,2)} + {l^2})|. (A.l) 
The expressions for twist-5 GDAs are identical up to the substitution 

f ^n7r~'~ pmz'^ t^ti-k'^ -^nTr"^ JC^'^^ A.^^~^ 'j^n-n'^ rpn7T~^ rpn7T~^ \ ^ 

\X '1 '2 '3 '2 '3 '2 '3 ^ T J 

r 07177^ rjnTT^ T/n7r^ T/n-TT^ /tn-n""'" /tn-n"^ rjiriTr^ rpnTr^ rrinTr^ T /' A 0\ 

1'52 1^2 1*^5 '^5 '^4 i-'s '-'4 '-'s J' l^-^J 

and the similar replacement of the neutron DAs on the r.h.s. 

Finally, the expressions for twist-6 GDAs are identical to twist-3 with the substitution 

and the similar substitution of the neutron DAs. 

2. ptt" distribution amplitudes 

Definition of the GDAs follows Eq. (|3.28p . For twist-4 we obtain 
Sr\l, 2, 3)— i{2Pf (1, 3, 2) - Sf{l, 2, 3)}, 

Pr\l, 2, 3)=-i{25f (1, 3, 2) - Pf (1, 2,3)}, 

Vr\l,2,3)^ iKf(l,2,3), Vr\l,'2,3) = ^Vi{l,2,3), 

Af (1, 2, 3)= i^^(l, 2, 3) , Af (1, 2, 3) = ^^^(1, 2, 3) , 

Tr"(l, 2, 3)=-iT|(l, 2, 3) , Tr\l, 2, 3) = ^ {r|(l, 2, 3) + 2rP(l, 2, 3)} , 

Tr\l, 2, 3)=i{rf (1, 2, 3) + 2r|(l, 2, 3)}. (A.4) 

The expressions for twist-5 and twist-6 GDAs are identical to those for twist-4 and twist-3, respectively, with 
obvious substitutions as described for the mr~^ case. 
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3. piT^ distribution amplitudes 



For completeness we present here the expressions for the p7r+ GDAs. They are defined as the corresponding 
matrix element of the trilocal Ught-ray operator 

Ol^l = e'^^v'^{a^zyp{a2z)u';{a^z) (A.5) 

and can be relevant, e.g. for weak form factors corresponding to the charged VF-boson exchange. Working 
out the necessary commutators we obtain 

[Qi,o^;^!;(z)]==-i{(75)aAO^;^!;(^) + {i,)pxdi%{z) + {i^UoiiHz)] , 

[Qldr^-{z)]^ '-{h5)axdtp-{z) + {i,)pxdl%{z) + {i,),xdlli{z)] , (A.6) 
In fact, a separate calculation is not needed in this case. For all GDAs: 

(1,2, 3)^-Sr* (1, 2, 3) + Sf{l, 2, 3) , 
(1,2, 3)=-Pr^ (1, 2, 3) + Pf (1, 2, 3) , 
(1, 2, 3)= (1, 2, 3) + F/(l, 2, 3) , 

Ar^(i,2,3)= Ar^(i,2,3)+Ani,2,3), 

Tf (1, 2, 3)^-7;""^ (1, 2, 3) + Tf (1, 2, 3) . (A.7) 



APPENDIX B: ^^-CORRECTIONS 



1. v^'^-"' corrections 



The method to calculate the 0{x^) contributions in the OPE of the baryon operators in a diquark-quark 
configuration (i.e. when the light-cone positions of two quarks coincide) is described in detail in Refs. (30l.[3^ 
and need not be repeated here. The functions V*^ can be found as solutions of the following equations for 
the moments: 



1 n 1 ,TTN,Al(d) , N 



1 



(n + l)(n + 2) 
1 

^ (n+ l)(n + 3) 
1 



(-2Vi + ^3 + ^4 + 21/5 



l(d)(n+2) 



(B.l) 



(n + l)(7i + 3) 
1 

^(n + l)(?i + 2) 



(-2^1 + 14 + V^4 + 2F5) 



(«)(n+2) 



(B.2) 



We present the solutions as 



V. 



■KN,M(d) 



(B.3) 
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where 



- X2f\^6 + 18x2 + 19x1 + lOfixl - 3(9 + lQfi)xl - 4/^(1 + 3x2 + 26x1 - 25a;^)|, 
= - .T2)^|6 + 18a;2 + U3xl + 297x1 ^ 408x| - 10A5'(1 - 3a;2)a;^ 

--2Vf{A + 12x2 + 7x1 + 599x1 - 612a;^)|, 
=-i|3x^(l - a;3)(87 - 33x3 - ISxj + llx^) + 156x^ ln[a;3] 

+ 2ff (1 -a;3)(4 + 4x3 -351a;^ + 169x^ + 490;;^- 55x1) -180x^ln[a;3] }, 
=-i|x^(l - X3)(1979 - 877x3 - 1513x^ + 2232x1 - 816x^) + 996x2 j^^j^^j 

- AV{^ (1 - X3)(4 + 4x3 + 1527x3 - 893x| - 923x;^ + 1613x^ - 612x^) + 720x3 ln[x3] |, 

= — ^(1 - X2)'''|2 + 6x2 - 75x^ + 43x^ + 4/i"(l + 3x2 + 26x^ - 25x^) 
v2 

+ 2/i*(4 + 12x2 + 89x^ - 55x^)|, 

= — ^(1 - X2)^|l0 + 30x2 + 425x^ - 337x^ - 30Vf (1 - 3x2)xl 
v2 

- 2A5'(4 + 12x2 + 17x^ + 569x^ - 612x^)|, 

= ^|(1 - X3)(-8 - 8x3 + 441x2 _ 239x^ - 59x;^ + 77x^) + 204x^ln[x3] 
v2 

+ 2(/i" + 2f^) (1 - X3)(4 + 4x3 - 351x^ + 169x^ + 49x;^ - 55x|) - 180x^ ln[x3] |, 

-- — ^{(1 - a;3)(16 + 16x3 + 171x^ - 941xi^ + 847x^ - 217x^) - 108x^ ln[x3] 
v2 

+ 4^5" (l-X3)(4 + 4x3 + 1527x^- 893x|- 923x^ + 1613x^-612x|) + 720x^ln[x3] |. (B.4) 



2. Aj^''^^ corrections 



Similarly wc get the equations for the A^^ functions: 

1 



jdx,xiAr^''^'\x,y 



(71 + l)(n + 2) 
1 

^ (7i + l)(n + 3) 



(-2^1+^3 + .44 + 2^5)'''^^"+'^ 



71- AT 



(B.5) 



(71 + 3)^1''^^"+') + (Ai - A2)('')("+2) - (Ai + A5)('^)("+i)' 



dX2 X2 Ai {X2) = 



1 



(71+ 1) (71 + 3) 
1 

"^(71+ l)(n + 2) 



{~2Ai +A3 + A4 + 2^)(")("+2) 



ttTV 



(B.6) 
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We obtain 

A =i^{^~ x^Yy^iD^ +fNDf j. (B.7) 

with the coefficients 

D^'^" =~]^^xl{29 - 45x2) + 2/i"(4 + 12x2 - 51x^ + Ahxl) - 4/{'(2 + 6x2 + 17x1 - 30x^)}> 
^u,pxO^_ 1 |g ^ 24a;2 _ 21x1 + 2^1 (4 + 12x2 - 8x2 - 609x^ + &l2xl) + Abxl + lQV^xl{l - 8x2)}, 
^«,«7r+^_l^|3^2('47 _ _ ^ju^2 + 6x2 + 17x^ - 30x^) + 2/f (4 + 12x2 - 221x^ + 255x^)|, 

^«,n7r+^_i_|-^g + 48x2 - 35x^ + 1323x^ - 1224x^ + 2Al{S + 24x2 - ISxj - 1233x^ + 1224x2) 

+ QV^{A + 12x2 - 9x^ - 609xi + 612x|)|. (B.8) 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 

[7] 

[8] 

[9] 
[10] 

[11] 

[12] 

[13] 
[14] 



N. M. KroU, M. A. Rudemian Phys. Rev. 93, 233 
(1954). 

Y. Nambu and E. Shrauner, Phys. Rev. 128, 862 
(1962). 

A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B 
36, 589 (1972). 

S. Scherer and J. H. Koch, Nucl. Phys. A 534, 461 
(1991). 

P. V. Pobylitsa, M. V. Polyakov and M. Strikman, 
Phys. Rev. Lett. 87 (2001) 022001. 
V. L. Chernyak and A. R. Zhitnitsky, JETP Lett. 
25, 510 (1977) Sov. J. NucL Phys. 31, 544 (1980); 
V. L. Chernyak, A. R. Zhitnitsky and V. G. Serbo, 
JETP Lett. 26, 594 (1977) Sov. J. Nucl. Phys. 31, 
552 (1980). 

A. V. Radyushkin, J INR report R2-10717 (1977), 
| arXiv: hep-ph/0410276 (English translation) ; 

A. V. Efremov and A. V. Radyushkin, Theor. Math. 
Phys. 42, 97 (1980) Phys. Lett. B 94, 245 (1980). 
G. P. Lepage and S. J. Brodsky, Phys. Lett. B 87, 
359 (1979); Phys. Rev. D 22, 2157 (1980). 

B. Pire and L. Szymanowski, Phys. Lett. B 622, 83 
(2005). 

J. P. Lansberg, B. Pire and L. Szymanowski, 
|arXiv:hep-p h/0610407 

N. Isgur and C. H. Llewellyn Smith, Phys. Rev. 
Lett. 52, 1080 (1984). 

N. Isgur and C. H. Llewellyn Smith, Nucl. Phys. B 
317, 526 (1989), Phys. Lett. B 217 (1989) 535. 
P. KroU, M. Schurmann and P. A. M. Guichon, 
Nucl. Phys. A 598, 435 (1996). 
A. V. Radyushkin, Phys. Rev. D 58, 114008 (1998). 



[is; 

[16] 
[17 

[is; 

[19 

[20; 

[21 

[22; 

[23; 

[24; 

[25; 

[26] 

[27 

[28' 
[2Q 

[3o; 

[31 



M. Diehl, T. Feldmann, R. Jakob and P. KroU, Eur. 
Phys. J. C 8, 409 (1999). 

K. Goeke, M. V. Polyakov and M. Vanderhaeghen, 

Prog. Part. Nucl. Phys. 47, 401 (2001). 

M. Diehl, Phys. Rept. 388, 41 (2003). 

A. V. Belitsky and A. V. Radyushkin, Phys. Rept. 

418, 1 (2005). 

A. V. Belitsky, X. d. Ji and F. Yuan, Phys. Rev. D 
69, 074014 (2004). 

M. Diehl, T. Feldmann, R. Jakob and P. KroU, Eur. 
Phys. J. C 39, 1 (2005). 

M. Guidal, M. V. Polyakov, A. V. Radyushkin and 

M. Vanderhaeghen, Phys. Rev. D 72, 054013 (2005). 

P. A. M. Guichon, L. Mosse and M. Vanderhaeghen, 

Phys. Rev. D 68, 034018 (2003). 

M. V. Polyakov and S. Stratmann, 

|arXiv:hep-p h/0609045 

I. I. Balitsky, V. M. Braun and A. V. Kolesnichenko, 
Nucl. Phys. B 312, 509 (1989). 
V. L. Chernyak and I. R. Zhitnitsky, Nucl. Phys. B 
345, 137 (1990). 

V. M. Braun, A. Khodjamirian and M. Maul, Phys. 
Rev. D 61, 073004 (2000). 

J. Bijnens and A. Khodjamirian, Eur. Phys. J. C 

26, 67 (2002).^ 

V. M. Braun, |arXiv:hep-ph/9801222| 

P. Colangelo and "K. Khodjamirian, 

larXiv:hep-p h/0010175 | 

V. M. Braun, A. Lenz, N. Mahnke and E. Stein, 
Phys. Rev. D 65, 074011 (2002). 
Z. G. Wang, S. L. Wan and W. M. Yang, Phys. Rev. 
D 73, 094011 (2006). 



21 



[32] Z. G. Wang, S. L. Wan and W. M. Yang, Eur. Phys. 

J. C 47, 375 (2006). 
[33] A. Lenz, M. Wittmann and E. Stein, Phys. Lett. B 

581, 199 (2004). 
[34] V. M. Braun, A. Lenz and M. Wittmann, Phys. Rev. 

D 73, 094019 (2006). 
[35] M. Q. Huang and D. W. Wang, Phys. Rev. D 69, 

094003 (2004). 
[36] M. Q. Huang and D. W. Wang, 

|arXiv:hep -ph/0608170, 

[37] Z. G. Wang, |arXiv:hep-ph/0609l55| 
[38] V. M. Braun, A. Lenz, G. Peters and 

A.V. Radyushkin, Phys. Rev. D 73 (2006) 034020. 
[39] J. D. Bjorken and S. D. Drell, Relativistic Quantum 



Fields (McGraw-Hill, New York, 1965). 
[40] V. Braun, R. J. Fries, N. Mahnke and E. Stein, Nucl. 
Phys. B 589 (2000) 381 [Erratum-ibid. B 607 (2001) 
433]. 

[41] B. L. loffe, Nucl. Phys. B 188, 317 (1981) [Erratum- 
ibid. B 191, 591 (1981)]. 

[42] V. L. Chernyak and I. R. Zhitnitsky, Nucl. Phys. B 
246, 52 (1984). 

[43] I. D. King and C. T. Sachrajda, Nucl. Phys. B 279, 
785 (1987). 

[44] V. L. Chernyak, A. A. Ogloblin and L R. Zhitnitsky, 

Z. Phys. C 42, 569 (1989). 
[45] J. Boiz and P. KroU, Z. Phys. A 356 (1996) 327. 



